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Universally Coupled Massive Gravity* * 
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Abstract 

We derive Einstein’s equations from a linear theory in flat space- 
time using free-field gauge invariance and universal coupling. The 
gravitational potential can be either covariant or contravariant and 
of almost any density weight. We adapt these results to yield uni¬ 
versally coupled massive variants of Einstein’s equations, yielding two 
one-parameter families of distinct theories with spin 2 and spin 0. The 
Freund-Maheshwari-Schonberg theory is therefore not the unique uni¬ 
versally coupled massive generalization of Einstein’s theory, although 
it is privileged in some respects. The theories we derive are a subset 
of those found by Ogievetsky and Polubarinov by other means. The 
question of positive energy, which continues to be discussed, might be 
addressed numerically in spherical symmetry. We briefly comment on 
the issue of causality with two observable metrics and the need for 
gauge freedom and address some criticisms by Padmanabhan of field 
derivations of Einstein-like equations along the way. 
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1 Introduction 

Constructing a relativistic gravitational theory based on principles such as 
an analogy to Maxwellian electromagnetism, the universal coupling of the 
gravitational field to a combined gravity-matter energy-momentum com¬ 
plex, and also the requirement that the gravitational field equations alone 
(without the matter equations, again in analogy with electromagnetic charge 
conservation) entail energy-momentum conservation was a major part of 
Einstein’s search for an adequate theory of gravity in 1913-1915 (see m,®)- 
Einstein subsequently downplayed these investigations [2j, and the above 
ideas later came to be associated with the non-Einsteinian field theory ap¬ 
proach to gravitation. Such a derivation of Einstein’s or similar gravita¬ 
tional field equations can use a priori preferred coordinates and a canonical 
energy-momentum tensor or a flat background metric and variational metric 
energy-momentum tensor, the difference being basically formal. 
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Several authors m-m have discussed the utility of a flat background 
metric r]^ u in general relativity or the possibility of deriving that theory, 
approximately or exactly, from a flat space-time theory^] A background 
metric allows introducing a gravitational energy-momentum tensor na , not 
merely a pseudotensor. As a result, gravitational energy and momentum 
are independent of the coordinates but dependent on the gauge [18]. If 
we want to regard the background metric seriously as a property of space- 
time and not just treat it as a useful fiction, then the relation between the 
effective curved metric’s null cone and that of the flat background must be 
considered. That issue was addressed with some success for the massless 
case of Einstein’s equations US- 

Preparatory to considering massive theories of gravity, we generalize our 
derivation of Einstein’s equations using gauge invariance for the free field 
and universal coupling [T2] to permit almost any density weight and either 
covariant or contravariant valency for the symmetric rank-two gravitational 
potential. This generality to some degree parallels that in Kraichnan’s clas¬ 
sic work [0], but our derivation has several improvements and can be easily 
adapted to massive theories. The choices of index position and density 
weight make no difference (after field redefinitions) in the massless theories, 
but they yield distinct massive theories here. Along the way, we address 
Padmanabhan’s recent objections to field derivations of Einstein’s equa¬ 
tions TV)]. 

Several authors recently discussed the subject of massive gravity with 
spin-2 and spin-0 components m , m- While it is permissible simply to 
postulate the nonlinear features (if any) of a mass term, it seems preferable 
to find well-motivated theoretical principles to constrain such choices. Some 
time ago, Ogievetsky and Polubarinov (OP) derived a two-parameter fam¬ 
ily of massive variants of Einstein’s equations |7j, which contains both our 
one-parameter families (quantization was considered briefly in [22]). Their 
derivation relied on gauge invariance (at least for the massless part of the 
Lagrangian density) but not universal coupling. Instead, they imposed a 
spin-limitation principle to exclude some degrees of freedom, many with 
the wrong sign and thus negative energy, from the full nonlinear interact- 

1 We note that general relativity was also derived from self-interaction on curved back¬ 
grounds m, m The possibility (or otherwise) of massless multigraviton theories anal¬ 
ogous to Yang-Mills theories was also investigated [16] . Moreover, it is interesting to 
consider (A)dS backgrounds, Lorentz-violating theories, and Chern-Simons topologically 
massive theories. But we consider only the most traditional form of the massive grav¬ 
ity problem here. We work in four space-time dimensions, but the derivations should 
generalize straightforwardly to any (integral) dimension not less than three. 
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ing theory. The Freund-Maheshwari-Schonberg (FMS) massive theory was 
originally derived using universal coupling with a canonical, not metric, 
energy-momentum tensor [23]. While the FMS theory is recovered among 
our results, it is not the unique universally coupled massive variant of Ein¬ 
stein’s equations, contrary to previous claims [23], [23]. Our derivation us¬ 
ing the metric energy-momentum tensor seems shorter and cleaner than the 
FMS derivation using the canonical energy-momentum tensor. In a future 
work, we will consider universal coupling using a tetrad, not metric, formal¬ 
ism, thus obtaining additional universally coupled theories, and will use the 
metric formalism to show that all the OP theories are universally coupled. 
All known universally coupled theories correspond to the OP family of the¬ 
ories; hence, their derivation and the metric universal coupling derivation 
currently lead to coinciding results. The FMS theory, subsequently adopted 
by Logunov and collaborators mi. also faces questions regarding positive 
energy and causality, on which we briefly comment in what follows. 

2 Effectively geometric theories from universal cou¬ 
pling and gauge invariance 

We previously derived Einstein’s theory and other effectively geometric the¬ 
ories using gauge-invariant free-field theories and universal coupling [ 12 ], 
which was significantly based on the work of Kraichnan [3] and Deser [9]. 
Whereas the gravitational potential was previously taken to be a symmet¬ 
ric rank-two covariant tensor field, we now generalize this derivation to the 
case of a density of almost any weight and either covariant or contravariant 
valency. Given the generality of Kraichnan’s derivation, it is not surprising 
that these generalizations again yield only Einstein’s and other effectively 
geometric theories. Nothing especially novel is obtained for massless theo¬ 
ries, but the derivations below are adapted to massive gravity for the first 
time. 

2.1 Free-field action for a covariant tensor density potential 

For the massless theories, an initial infinitesimal invariance (up to a bound¬ 
ary term) of the free gravitational action is assumed. For the subsequent 
derivation of massive theories, the gauge freedom is broken by a natural 
mass term algebraic in the fields, but the derivative terms retain the gauge 
invariance. 

Let S{ be the action for a free symmetric tensor density 7 ^ (of density 
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weight —l, where l 7 ^ 1 / 2 ) in a space-time with a flat metric tensor q^ u in ar¬ 
bitrary coordinates. The torsion-free metric-compatible covariant derivative 
is denoted by d hence, d a q lw = 0. It is convenient to use not the flat metric 
itself but its related densitized metric fj^ = q of weight —l. We 
note that in the forbidden case l = 1 / 2 , q^ v is noninvertible: q^u{\/—q ) -1 / 2 
determines only the null cone, not a full metric tensorJl The field 77 ^ turns 
out to be the gravitational potential. Although it has become customary 
in work on loop quantum gravity to denote the density weight of fields by 
placing the corresponding number of tildes above or below a symbol to ex¬ 
press its positive or negative density weight, that custom is impossible here. 
The density weight l (Z 7 ^ 1/2) can be large, nonintegral, or even irrational, 
and we therefore merely write a tilde over most densities. All indices are 
respectively raised and lowered with q^ v and q fiu with two exceptions. For 
the densitized flat metric q^ v , the oppositely densitized inverse flat metric 
is 77 ^". Similarly, the inverse of the densitized curved metric g defined be¬ 
low is g ^ v . Because we use tensor densities extensively, we recall the forms 
of their covariant and Lie derivatives. A (1, l)-density qbg of weight w is a 
representative example. The Lie derivative is given by [25] 

£* 4 % = +<^0 +«'$£%, ( i ) 

and the r/-covariant derivative is given by 

( 2 ) 

where T^ are the Christoffel symbols for q /1U . After the curved metric g^ w 
is defined below, the analogous ( 7 -covariant derivative V with the Christoffel 
symbols follows. 

The desire to avoid ghosts motivates gauge invariance for linear theo¬ 
ries [26] . We require that the free field action S{ change only by a boundary 
term under the infinitesimal gauge transformation 

'J/j.v * qin- 1 T 1 >,vs fi'Y/j.v — , + du T cq^i/d //>, (3) 

where c 7 ^ — 1/2 and is an arbitrary covector density field of weight — zJl 
We can expect the appearance of a connection between Z and c. For any 
S'f invariant in this sense under @, a certain linear combination of the free 

2 This exceptional case was of interest in deriving slightly bimetric theories m where 
yj—g can appear in the field equations of the interacting theories. 

3 The case c = —1/2 merely gives a scalar theory in the somewhat comparable work of 

OP [ 3 . 
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field equations is identically divergenceless, as we now show. The action 
changes by 


5S { 



' 5Sf 
_S%i> 


(di + d + crj^yd 01 ^) + 


/ 


d 4 5 x /%. 


(4) 


The explicit forms of the boundary terms given by and are not 
needed for our purposes. Integrating by parts, letting have compact 
support to annihilate the boundary terms, and using the arbitrariness of 
we obtain the identity 


d„ 


5S{ c av SSf 




(IV 


8% 


= 0. 


(5) 


This is the generalized Bianchi identity for the free theory, which in the most 
common case is a linearized version of the original geometric Bianchi identity. 
A natural choice for S'f is the linearized GR Lagrangian density, but no such 
detailed assumptions on the form of S'f are used here. Because noninteracting 
sourceless fields are unobservable [27], the theory is interesting only after an 
interaction is introduced. 

As Kretschmann pointed out long ago in response to Einstein, any theory 
can be given a generally covariant formulation in the sense that the equations 
hold in any coordinate system, Cartesian or otherwise m-m- The result¬ 
ing formulation might be called weakly or trivially general covariant. Often 
achieving weak general covariance involves using absolute objects m, [so], 
such as a flat metric tensor, as it does here. The distinctively novel aspect 
of Einstein’s theory of gravity is supposedly its lack of absolute objects or 
“prior geometry” [29], [32] FI This property might be called strong or non¬ 
trivial general covariance. This distinction between two senses of general 
covariance was discussed previously [35]. The derivation presented here and 
in some other works starts with a weakly generally covariant theory and 
concludes with a strongly generally covariant one without absolute objects H 

4 Actually, the matter is more complicated: R. Geroch and Giulini recently noted in 
effect that g, the determinant of the metric tensor, is an absolute object because any point 
has a neighborhood with a coordinate system such that the component of g has the value 
— 1 [33], [34]. 

5 The fact that g counts as absolute in general relativity in the Anderson-Friedman 
absolute-object program suggests that general relativity is not strongly covariant after all. 
But the point remains that two very different notions of general covariance are in play. In 
addition to the field equations, the topology, boundary conditions, and causality should 
be examined in seeking absolute objects m We emphasize that the massless cases are 
considered here; obviously, the massive theories contain the flat metric and are therefore 
not strongly generally covariant. 


5 








Therefore, the quite misleading claim m that flat space-time derivations 
of Einstein’s equations only result in general covariance because they feed 
it in at the beginning should not be accepted. The two senses of general 
covariance are very different; Padmanabhan’s criticism commits the fallacy 
of equivocation. 


2.2 Metric energy—momentum tensor density 

If the energy-momentum tensor is to be the source of the gravitational 
potential 7 ^, then consistency requires that the total energy-momentum 
tensor be used including the gravitational energy-momentum, not merely 
the nongravitational (“matter”) energy-momentum, because only the total 
energy-momentum tensor is divergenceless in the sense of d u [9] or, equiva¬ 
lently, in the sense of a Cartesian coordinate divergence. To obtain a global 
conservation law, a vanishing coordinate divergence for the four-current is 
needed. 

An expression for the total energy-momentum tensor density can be 
derived from S using the metric recipe s, m, m, [as] as follows. The 
action depends on the flat metric density fj /JjV , gravitational potential 7 ^, 
and matter fields u. Here, u represents an arbitrary collection of dynamical 
tensor fields of arbitrary rank, index position, and density weight. Using the 
OP spinor formalism that uses the “square root of the metric” rather than 
a tetrad or other additional structure mi we can likely also include spinor 
fields. 

Under an arbitrary infinitesimal coordinate transformation described by 
a vector field the action changes by the amount 


[ a [5S 6S 5S \ , . 

5S = J d A x £^u + -rr— £$ri tiv + 3 % J , (6) 

with boundary terms from g vanishing because is compactly sup¬ 
ported. But S' is a scalar, and hence 5S = 0. Letting the matter and 
gravitational field equations hold, integrating by parts, discarding vanishing 
boundary terms, and using the arbitrariness of the vector field gives 


d u 



:Vapn 


JUS 



= 0 . 


(7) 


This quantity is an energy-momentum tensor density for matter and gravita¬ 
tional fields. The flatness of r]^ v is relaxed in taking the functional derivative 
SS/Srjftu and is restored later. This move is sometimes criticized [T9] , but it 
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is unobjectionable even in flat space-time because it is merely a formal trick 
useful for defining the energy-momentum tensor, not an illicit use of curved 
space-time. Using the connection between the Rosenfeld energy-momentum 
tensor and the symmetrized Belinfante canonical energy-momentum ten¬ 
sor [136], [3H], we could regard the metric recipe as a mathematical shortcut 
to the conceptually unimpeachable but mathematically inconvenient Belin¬ 
fante tensor modified with terms proportional to the equations of motion. 
The metric energy—momentum tensor is not unique at this stage, because 
terms proportional to the equations of motion and their derivatives, as well 
as superpotentials, can be added. Another option would be to use Lagrange 
multipliers and let the background metric be flat only on-shell [39j. As is 
shown below, using the superpotential freedom judiciously is important in 
deriving the field equations. 


2.3 Full universally coupled action 


We find an action S satisfying the plausible physical postulate that invert¬ 
ible linear combinations of the Euler-Lagrange equations are just invertible 
linear combinations of the free field equations for Sf augmented by the total 
energy-momentum tensor. A simple way to impose this requirement reduces 
to the condition 

SS _ SSf SS 

Sl/iv S% v 8fj,j, v ’ 

where A = — \Z32ttG. (Because the linear combinations of the free Euler- 
Lagrange equations and of the full Lagrange equations with interaction are 
invertible in this case, the linear combination can now be applied to the 
energy-momentum tensor.) It seems prudent to set c = —l to make the 
generalized Bianchi identity and the energy-momentum tensor density take 
similar forms. 

The basic variables here are the gravitational potential 77^ and the flat 
metric density fj^ u . But we can freely change the variables in S from 77^ 
and fj^ to the bimetric variables g^ w and fj^ [I], where 


9fiv — V/w (9) 

(We can then define the metric g^ u from g^ v using matrix algebra and then 
define the (7-covariant derivative V as usual, but we have little need to use 
V explicitly.) 

Equating coefficients of the variations gives 

SS _ SS 

Sfjfn/ Sfj^u 


9 + 


SS 


S~9 


LIV 


( 10 ) 


7 











for drj^ and 


( 11 ) 


6S _ 6S 

filin' ^9[IV 

for d'y^. Combining these two results gives 


A 


SS 

Sfj^ 


7 = A 


8S 


H 


'fils 


5S 

/Mis’ 


( 12 ) 


Equation (1121) splits the energy-momentum tensor into two parts: one that 
vanishes when gravity is on-shell and one that does not. Using this result 
in universal-coupling postulate (JH) gives 


A 


8S 


SS { 



(13) 


Taking the divergence, recalling free-theory Bianchi identity ([5|h and 
using c = —l, we derive 


d„ 


5S 


5fj, 


11V 


l 


5S 


~ 9 ~2 n ^ 


5 =0. 


(14) 


The quantity in parentheses is exactly {yf—rj) 1 SS/Srj^g. Hence, the part 
of the energy-momentum tensor not proportional to the gravitational field 
equations has identically vanishing divergence (on either index), i.e., is a 
(symmetric) “curl” [29]. The splitting of the energy-momentum tensor 
ensures that in the massless case, the gravitational held equations alone , 
without separately postulating the matter equations, entail conservation of 
energy-momentum for the resulting effectively geometric held equations. 

Because the quantity SS/Srj^lg is symmetric and has identically vanish¬ 
ing divergence on either index, it necessarily has the form [50] 


9 = ^d p d a (M [pp][rru] + M [l/p][ap] ) + (15) 

where is a tensor density of unit weight and B is a constant. This 

result follows from the converse of the Poincare lemma in Minkowski space- 
time. We cannot choose JV[P pav arbitrarily but must choose it such that the 
term 5Sf/5 ; y f j jl/ is accommodated. The freedom to add an arbitrary curl must 
therefore be used in a quite definite way. As Huggins, a student of Feynman, 
showed in his dissertation m and Padmanabhan recently emphasized m, 
coupling a spin-2 held to the energy-momentum tensor does not lead to a 
unique theory, because of terms of this curl form. Rather, as Huggins said 


8S 

8r]uu 
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(p. 39 in [31]): “an additional restriction is necessary. For Feynman this 
restriction was that the equations of motion be obtained from an action 
principle; Einstein required that the gravitational field have a geometric 
interpretation. Feynman showed these two restrictions to be equivalent.” 

Gathering all dependence on rj^ (with g^ independent) into one term 
yields S = S\[g^ v ,u] + S 2 [g fl u, 77^, u]. It is easy to verify that if 


5 2 = ^ J d 4 x Repair]) MV’ vp(T (r/^, < 7 ^, u) + J d 4 xa%+ 2 H J ^xy/^rj, 

(16) 

then 5 S 2 /Sg^lg has exactly the desired form, while S 2 does not affect the 
Euler-Lagrange equations, because 882 /begin = 0 and 5 S 2 /SU = 0 identi¬ 
cally [3]0 The coefficient B of the four-volume term is naturally chosen to 
cancel any zeroth-order term (such as from a cosmological constant) in the 
action such that the action vanishes when there is no gravitational field. 
The four-divergence resolves worries [19) about obtaining terms that 
are not analytic in the coupling constant A. It is unclear whether the Hilbert 
action is best in any event, given its badly behaved conservation laws [32] , 
Hence, the universally coupled action for the covariant tensor density 
case is 


(17) 

The boundary term is at our disposal; if is a unit-weight vector density, 
then 5 is a coordinate scalar. Using the effective curved metric density 
g^u, we can define an effective curved metric by g jlu = g^{y/—g )and an 
inverse curved metric density g^ v = g^ u {y/—g ) l ■ 

For 5i, we choose the Hilbert action for general relativity plus minimally 
coupled matter and a cosmological constant: 


Si = 


1 


167rG 


J d A x ^gR(g) - J d A x + 5 ma tt er [g,iv , u ]. (18) 


It is well known that the Hilbert action is the simplest (scalar) action that 
can be constructed using only the metric tensor. If the gravitational field 
vanishes everywhere, then the gravitational action should also vanish. In the 

6 If it seems that using this term (1/2) f d A x R^ vpa {ri)M >lvp! ’gu) is too clever 
to be invented without knowing Einstein’s theory in advance and thus cheating [19], then 
a symmetric curl term d p d a (M^ ,lp ^ al '^ + Ai^ p ^ <T ^)/2 + B^/—gg^ should be simply added 
to the metric energy-momentum tensor by hand using the usual underdetermination of 
the energy-momentum tensor. 
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massless case, the result is that B = A/(167 tG). For the generalization to 
the massive case considered below, the gauge-breaking part of the mass term 


introduces another zeroth-order contribution that also needs to be canceled. 


It is also possible to couple the matter to the Riemann tensor for g^ u or 
to allow higher powers of the Riemann tensor in the gravitational action, if 
desired. In the massless case, we could set A = 0 [23] . 

3 Massive universally coupled gravity for a covari¬ 
ant tensor density potential 

Our goal is to generalize the above derivation to obtain one or more massive 
finite-range variants of Einstein’s equations. Such field equations would be 
related to Einstein’s much as Proca’s massive electromagnetic field equa¬ 
tions are related to Maxwell’s. But a spin-2-spin-0 massive theory would 
have ghosts at the level of the free linear theory. Good linear behavior is 
generally required as a guide to good nonlinear behavior. But good linear 
behavior seems neither necessary nor sufficient for good nonlinear behav¬ 
ior. In the present case, it seems quite possible that the nonlinear form of 
the Hamiltonian constraint cures the bad behavior of the linear theory. We 
briefly discuss this matter below. 

It can be expected that the mass term for a free field is quadratic in 
the potential and lacks derivatives. The free-held action Sf is now assumed 
to have two parts: a (mostly kinetic) part Sjq that is invariant under the 
previous gauge transformations as in the massless case above and an alge¬ 
braic mass term Sf m that is quadratic and breaks the gauge symmetry. We 
seek a full universally coupled theory with an action S that has two corre¬ 
sponding parts: S = So + S ms . They are the familiar part So (yielding the 
Einstein tensor, the matter action, a cosmological constant, and a zeroth- 
order four-volume term) and the new gauge-breaking part S'ms, which also 
has another zeroth-order four-volume term. As it turns out, the mass term 
is constructed from both the algebraic part of So (the cosmological constant 
and four-volume term) and the purely algebraic term S ms . 

Requiring Sfo to change only by boundary terms under the variation 
liiv —> ly.v + d^v + dvin + cr hj.vd a £a for c/-1/2 implies the identity 



( 19 ) 
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We again postulate the universal coupling in the form 


5S _ 5S{ 5S 

Slixv 5%u Srjn V ' 


( 20 ) 


Changing to the bimetric variables g^ v and fj^, as before, implies that 


5S{ 


5S 




fiv 


( 21 ) 


We now introduce the quantities Sf m and 5 ms , Sf = Sfo + Sfm and S = 
So + 5ms, in order to treat the pieces that existed in the massless case 
separately from the new pieces in the massive case. We thus obtain 


SSio 


SSfrn 


5So_ 

SfjflU 


g + X 


^ms 

5rj^ 


a- 


( 22 ) 


Assuming that the new terms Sf m and S ms correspond, we separate this 
equation into the familiar part SSfo/Sj^ = X(5So/5f] flu )\g and the new 
part SSf rn /S^'^ = X (5S ms /5fj^ u )\g. Using invariance (fl9l) . we derive the 
form of So as 

So = Slight,, u] + S 2 , (23) 


as in the massless case. We again choose the simplest case and obtain the 
Hilbert action with a cosmological constant, with matter coupled only to 
the curved metric. 

The new part in the massive case is 


£S fm 


, 5S ms 


Assuming that the free-held mass term is quadratic in the gravitational 
potential, we find that its variational derivative is 


5Sfo 


= a y /^rj ; y a p(f) o ‘ ,1 rj 01 ' + bfj 010 fj^). 


fllS 


Changing to the bimetric variables gives 
a V~9, 


X 


<-g a g + rjapW^f, 01 ' + bf, a0 f,n = A 


SSry 






9- 


We take the expression for S ms in the natural form 

S ms = j d A x ( pg a/ 3 fj a0 + q)V^, 


(24) 


(25) 
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where p and q are real numbers to be determined below. We note that the 
term \J—g itself, which gives a cosmological constant, plays no role here and 
is already included in So- Using the relation \J—r\ = (■>/—fj) 1 ^ 1 ~ 2l \ we cal¬ 
culate (SS ms /Sfi llu )\g. Equating A (SS ms /5fi, lu )\g to i.e., equating 

corresponding coefficients, determines several of the constants. Equating the 
coefficients of the y/—rjfj terms gives q = (2 — 4Z)a(l + 46)/A 2 , equating 
the coefficients of the \J—rj ff' u fj a ' g a p terms gives p = —ab(2 — 41)/X 2 , and 
equating the coefficients of the 'J—r}ff lil ff v g a p terms gives p = a/ A 2 . Using 
the last two results together gives b = 1/(41 — 2). Using all three results 
together gives q = —2a(2/ + l)/A 2 . Combining the algebraic piece of So with 
S ms gives 

S a i g = ~-^Q J d A xy^g+2B J d 4 x^fj+-^ j d 4 x (g afi fj a0 -41-2)^/^. 

(26) 

When the gravitational potential vanishes, S a i g and hence the zeroth- 
order term should also vanish. Imposing this condition gives 

A a(l — 21) 

167 tG A 2 

Because our goal is to find a massive generalization of Einstein’s theory, 
not a theory with an effective cosmological constant, we require that the 
first-order term in 7 ^ also vanish. Because A 2 = 32irG, it follows that 
A = a(l — 21)/2. Hence, the sign of the formal cosmological-constant term 
depends on the density weight of the initially chosen potential. We also 
expect the quadratic part of the algebraic component 5 a i g of the action to 
agree with the free-field mass term Sf m . After a binomial expansion and 
some algebra, we see that this is the caseQ The weak-field expansion of the 
full massive nonlinear action S allows relating the coefficient a to the mass 
m of the spin-2 gravitons: o = —m 2 . For nontachyonic theories, we impose 
the condition a < 0 . 

Combining all these results gives the total massive action S, which de¬ 
pends on the spin -2 graviton mass m and the parameter l controlling the 

7 In [24], there is a mistake in the binomial expansion for y/—g between Eqs. (43) 
and (44). 
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relative mass of the spin -0 ghost: 


S 


16+ ‘S'matter \dfiv i 'A T 

+ \J d 4 xR^ upa (ri)M fJ ' l ' p,T [f} llu ,g^, u\ + j d 4 xa M , M - 


9 

m 

I 67 tG 


d 4 x 


V-a( 2 l - !) - V~~y( 2 i +1) + 


( 27 ) 


where l A 1/2- These theories are all universally coupled, contrary to the 
claim that only the FMS theory has this property [23], [24j . 

The Euler-Lagrange equations are easily found if the metric g fJV is used 
as the dynamical variable. The result is 


5S 


V=9 C nv m2 
167tG 167tG 


21 - 1 


V^99^+ 


(V=v ) l+1 

4 (>/=?)' 


{2rT 


lTf*9c0Sn 


(28) 

Following Boulware and Deser [24] . we can linearize these theories and 
check whether the spin-0 component is tachyonic. Nontachyonicity entails 
— 1/2 < l < 1/2. For l = 0, the spin-0 ghost has the same mass as the 
spin -2 degrees of freedom, and this theory (with density weight zero for the 
potential) is hence the cleanest in the set. (The connection between density 
“weight” and graviton “mass” is an amusing linguistic accident.) Investi¬ 
gating various masses for the spin -0 degree of freedom might have some 
empirical consequences. This is especially the case in large-scale and homo¬ 
geneous situations, for example, in cosmology [21j. The ratio of the spin-0 
mass mo to the spin -2 mass m 2 is given by 


+ 


d Smarter 


?71g _ — 4 1 2 + 1 
m,n 2 1 2 + 1 


(29) 


Hence, the mass is an even function of l. For l = —1/2, the spin-0 degree of 
freedom is massless, and \J—r) is absent from the mass term. As l —> 1/2, 
the scalar again becomes light, and the coefficient of y / —g tends to zero, 
although the value l = 1/2 is forbidden. Between these massless endpoints, 
the spin-0 degree of freedom becomes heavier. At the midpoint l = 0, 
the scalar has the same mass as the spin -2 field, giving a simple form of 
the wave equation for the linearized massive Einstein equations. Thus, the 
spin -0 ghost is never heavier than the spin -2 degrees of freedom, and the 
1^0 theories hence have weaker gravitational attraction at large distances 
and in homogeneous situations. 
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4 Derivation of a contravariant tensor density po¬ 
tential: Massless case 


We briefly present the contravariant analogue of the above derivation of Ein¬ 
stein’s equations. The gravitational potential is now a contravariant sym¬ 
metric tensor density field 7 ^" of density weight l, where l 7 ^ 1/2. In addition 
to the obvious moves of some indices, we introduce some sign changes in the 
contravariant case. The hope for a simple rule relating index moves and sign 
changes is disappointed because the symmetry between the ( 0 , 2 ) theory of 
weight —l and the (2,0) theory of weight l is broken: the Lagrangian den¬ 
sity is a scalar density of weight 1, not weight 0. The consequences of this 
imperfect symmetry are more apparent in the massive case below than in 
the current massless case. It is convenient to use not the inverse flat metric 
itself but the densitized related metric fp iV = rf u (^J—r) f of weight l, where 
l / 1 / 2 . 

In the massless theories, we assume an initial infinitesimal invariance (up 
to a boundary term) of the free gravitational action Sf under the infinitesimal 
gauge transformation 7 ^ —► 7 ^ + S'y ^ 1 ', where 

Vf' = + d u e + crf v dj a , (30) 


c 7 ^ —1/2, and is an arbitrary vector density of weight l. It proves 
expedient to set c = — l. For any Sf invariant in this sense, a certain linear 
combination of the free-held equations is identically divergenceless: 

[)!<■ ( _ L r , ^ — 0 (31) 

d 5r a )~ 1 j 

This is the generalized Bianchi identity for the free-held theory. Local 
energy-momentum conservation, which holds with the use of the Euler- 
Lagrange equations for the gravity 7 ^" and matter u, can be written as 


d u 



l ~ad ~ 
Vuv 



= 0 . 


(32) 


We write the universal-coupling postulate in the form 

6 S _ _6 S l 5S 

firy/U' () lj^ l/ 


(33) 


(the reason for choosing of the sign of the term containing the energy- 
momentum tensor soon becomes clear). We obtain S, changing from 7 ^" 
and fj to the bimetric variables g 111 ' and , where 


= ff u + 


(34) 
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The coefficient of A 7 is chosen such that the covariant and contravariant 
cases, as far as is easily achieved, define the gravitational potential similarly. 
In particular, to the linear order in the potential in Cartesian coordinates, 
the traceless part of the gravitational potential 7 has the same observable 
significance, whether it is the density difference between the covariant curved 
and flat and curved metrics or the density difference between the contravari¬ 
ant flat and curved metricsll Equating coefficients of the variations gives 


5S 

5fj^ 


5S 


7 = 


5f]f iv 


9 + 


5S 

5g^ u 


(35) 


and 


whence we obtain 


6S _ x 6S 
~ 5g^' 


5S { 


-A 


5S 


5fj^ v 


(36) 

(37) 


The use of the generalized Bianchi identity implies that S splits into a com¬ 
ponent it] and a component S 2 that takes form (fl6l) . The quantity 

1 S 2 contains all the ineliminable dependence on the background metric and 
does not contribute to the field equations. The simplest choice of S± gives 
the Hilbert action for Einstein’s equations and a cosmological constant, as 
in the covariant case. The specific choice from the allowed values of l makes 
no difference in the massless case. 


5 Derivation for a contravariant tensor density po¬ 
tential: Massive case 

The choice of the density weight l makes a difference in the massive gen¬ 
eralization of this contravariant derivation, much as in the covariant case. 
The FMS-Logunov theory turns out to be the 1=1 contravariant universally 
coupled massive theory. While in some clear senses the 1=1 theory is the 
best of the contravariant massive theories, it is not the only such theory. 
The free-field action Sf again has two parts: the part Sfo that is mostly 

8 If = 77 v — X'jfi.v, then the inverse metric yields an infinite series expansion (at 
least formally), whose first term has a sign different from what a naive index raising might 
suggest: < 7 „ = r/^ + + ■ ■ • implies that Ri Ay^", not If the exact 

relation g 11 " — tj 111 ' = Xi holds (the sign of the A term is important), then 1/7 v ~ 7 M i/. 
Thus, the meaning of the gravitational potential is insensitive to a sign change, and it is 
therefore easier to compare the various massive theories. 
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kinetic and has a local gauge symmetry and an algebraic mass term 5f m . 
The action S of the full theory again splits into two parts, Sq and S ms . 

Requiring Sfo to change only by a boundary term under the infinitesimal 
variation 

5^ v = + 8 v e - IrT d a i a , (38) 

where l ^ 1/2, implies the identity 

dfl f _ 1 \ = 0 ( 39 ) 

\8Y V 2 <57 aOL ) 

We again postulate the universal coupling in form (l33l) . We change to the 
bimetric variables g^ v and fj fiU . Letting the new mass terms and the terms 
previously present agree separately, we find that the mass terms satisfy 


&Sfm 


(40) 


The action Sf m is chosen to be quadratic in the gravitational potential and 
to satisfy 

SSim = aV-Vl^i.VanVpv + bfj a pf}^). (41) 


S'y fU ' 


The quantity S ms is naturally chosen in the form 

S ms = I d 4 x {pg 01 ^fja/3 + q)'/~V 
for unspecified p and q. Matching the coefficients of several terms gives 


(42) 


q = 


b = — 


— (2 — 4£)a(l + 46) 


A 2 


1 


4 / — 2 ’ 


P = 


q = 


ab(2 - 4 1) 
A 2 ’ 

2a(2l - 3) 
A 2 ' 




Requiring the zeroth-order algebraic term in S to vanish gives 

p A q(i — 2/) 

16v tG A 2 ' 


Requiring the first-order algebraic term to vanish, after some algebra, gives 
A = —a(l — 21)/2. The second-order term agrees with the free-held mass 
term, as could be hoped. The spin-2 graviton mass m is given by a = —m 2 . 
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Combining all these results, we obtain the total massive action S, which 
depends on the spin-2 graviton mass m and the parameter l controlling the 
relative mass of the spin-0 graviton: 


S = 


1 


I 67 tG 


J d 4 x \f—g R(g) + S matter [g^,u\ + 


+ \j d 4 x R, W!>a (ri)M IJl ' pa [rf JV , , u] + j d 4 x — 


m 


16irG 


d 4 x -s/^g(2l - 1) + v / 3 ??(2/ - 3) + g^fj^ 


(43) 


where l ^ 1/2. The empirically doubtful Fierz-Pauli mass term is not among 
the universally coupled theories considered in this paper. All these theories 
are contained in the OP 2-parameter family. 

While we could use some contravariant tensor (or tensor density) to find 
the Euler-Lagrange equations, the equations are more easily compared with 
those previously found if the metric is used. The resulting equations are 


SS 


\Z~~9 m2 

167tG 167tG 


1 - 21 
2 


^9 9^ W + 


(V^9) l Vaf3 

4 (v^y- 1 





(44) 

Linearizing these theories shows that the spin-0 field is not a tachyon if and 
only if 1/2 < l < 3/2. For 1 = 1, the spin-0 ghost has the same mass as 
the spin-2 degrees of freedom, and these field equations, which correspond 
to the FMS-Logunov theory, are hence the cleanest in the set. At the linear 
level, the (0, 2) theory of weight —l is identical to the (2, 0) theory of weight 
l + 1. Hence, for the (2,0) theory of weight l + 1, the ratio of the spin-0 
mass m.o to the spin-2 mass 777-2 is given by 


+ 


^matter 

^9^u 


777-g —4(Z — l) 2 + 1 

777-2 2(Z — l) 2 + 1 


The mass is an even function of l — 1. For l = 3/2, the spin-0 degree 
of freedom is massless, and y/—r] is absent from the mass term. As l —> 
1/2, the scalar again becomes light, and the coefficient of \J—g tends to 0, 
although the value l = 1/2 is forbidden. Between these massless endpoints, 
the spin-0 degree of freedom becomes heavier. At the midpoint l = 1 (the 
FMS-Logunov theory), the spin-0 graviton has the same mass as the spin-2 
graviton, giving a simple form of the wave equation for the linearized massive 
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Einstein equations. Hence, the spin-0 ghost is never heavier than the spin-2 
degrees of freedom for this family of theories, although such can occur for 
the larger family of OP massive theories [7]@ 

6 Massive gravities and experiment 

In both experimental [43) and theoretical contexts, it is common to speak of 
the mass of the graviton, as if all gravitons must have the same mass. While 
all gravitons do have the same mass in the most famous spin-2-spin-0 mas¬ 
sive gravity (developed by FMS and studied by Logunov and collaborators), 
the existence of the OP theories shows that massive gravity has two mass 
parameters that should be tested experimentally. It would be worthwhile 
to ascertain to what degree the tacit assumption of equal spin-2 and spin-0 
masses is actually used in finding experimental bounds on massive gravity. 
Astrophysical tests for changes in the behavior of the degrees of freedom 
present in massless general relativity primarily bound the spin-2 mass. The 
empirical bounds on the spin-2 graviton mass are so tight that the spin-2 
part of the mass term is empirically negligible except in strong fields or over 
cosmic distances m, ED; these regimes are also those investigated by Lo¬ 
gunov and collaborators. The flexibility in the spin-0 mass in the massive 
theories derived here opens some phenomenological opportunities by increas¬ 
ing the range of the spin-0 repulsion that counterbalances some of the spin-2 
attraction. The larger family of OP massive theories permits either longer 
or shorter range for the spin-0 repulsion compared with the spin-2 attrac¬ 
tion. Babak and Grishchuk recently investigated a similar phenomenological 
flexibility [21] . They noted that their massive spin-2, massless spin-0 spe¬ 
cial case agreed with general relativity in cosmological contexts because of 
the high degree of symmetry. Hence, cosmological limits on the graviton 
mass(es) primarily bound the spin-0 mass. While the OP theories have non¬ 
linear mass terms motivated from first principles in contrast to Babak and 
Grishchuk’s linear mass terms motivated by mathematical simplicity, simi¬ 
lar qualitative behaviors of the two kinds of two-parameter massive gravities 
can be expected outside highly nonlinear regimes. 

In view of the tight empirical bounds on the graviton masses, the observ¬ 
able consequences of a mass term are rather difficult to detect. But there 
are two important theoretical issues that arise at the classical level for mas¬ 
sive variants of Einstein’s equations. The first is the well-known question of 

9 It is reassuring that in the conformally flat special case g the (2, 0) theory 
of weight l + 1 and (0, 2) theory of weight —l coincide. 
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stability, positive energy, etc., given the wrong-sign spin 0. The second is 
the question of causality: massive gravity is a special relativistic field theory 
with rjfnj observable, but there is reason to fear that field propagation might 
violate causality by having the light cone of the effective metric leak 
outside that of rj^. We now turn to these issues. 

7 Positive energy 

It has long been argued that massive variants of Einstein’s equations pose 
the unpleasant dilemma [24) . [44] that either the mass term is of the Fierz- 
Pauli form with 5oo 3 degrees of freedom (pure spin 2) and is empirically 
falsified by having a discontinuous massless limit or the theory has 6oo 3 
degrees of freedom including a wrong-sign spin 0 (a spatial scalar density) 
and instability arises after linearization. The former problem is the van 
Dam-Veltman-Zakharov discontinuity [35]) about which a large literature 
has appeared in the last few years after a long period of relative quiet. More 
relevant for our purposes is whether the massive theories with 6oo 3 degrees 
of freedom (spins 2 and 0) are unstable. 

Contrary to widely held views, Visser argued that the massive theories 
with 6oo 3 degrees of freedom might well be stable m- More recently, Babak 
and Grishchuk argued that such theories actually are stable (21] . Concerning 
the specific case of the FMS theory, the authors of the theory were them¬ 
selves unconvinced of instability [23], although they did not follow up on the 
matter after arguments for instability were published. In the middle to late 
1980s, Logunov and collaborators (such as Loskutov and Chugreev) adopted 
the FMS theory as the massive version of the relativistic theory of gravity. 
They argued that this specific theory might well be stable im, lineariza¬ 
tion arguments notwithstanding. More compellingly, Loskutov calculated 
the gravitational radiation from a bounded source and concluded that it is 
in fact positive definite [36], even though the theory has a wrong-sign spin-0 
component. It is curious that this conclusion has received so little response. 

While the question of positive energy (or positive mass, as is often said 
in a gravitational context) has not been settled with a favorable outcome (in 
the sense of a general proof that all exact solutions satisfying certain energy 
conditions have positive mass), neither do the arguments for instability from 
linearization seem compelling. It is useful to show using a Hamiltonian 
formalism that linearization is untrustworthy because it essentially changes 
the form of the Hamiltonian constraint such that instability becomes more 
plausible than is true for the exact nonlinear theory. Although the wrong- 
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sign spin 0 is present in the nonlinear theory, its kinetic energy is related 
to that of the right-sign degrees of freedom such that it can easily radiate 
only together with the positive energy degrees of freedom. Precisely this 
feature is lost upon linearization. This important feature of the Hamiltonian 
constraint depends essentially on a cubic term in the Hamiltonian density 
and hence on a quadratic term in the field equations. 

In the case of the FMS-Logunov contravariant weight-1 theory, the field 
equations impose a lower bound on the Hamiltonian density a little below 
zero, in contrast to boundary terms (which should be annihilated, at least in 
static contexts, because of the exponential Yukawa falloff of the gravitational 
potential of bounded systems). Discarding unimportant terms and factors 
from action (1431) . we obtain the Lagrangian density 

£ = V^gR(g) ~ rn 2 (-^g- y/=T) + ^V^gg lxu ri f ,\ (45) 

We use the ADM (3+l)-dimensional splilP^I |32] and choose coordinates 
(Cartesian, spherical, or the like) such that 7700 = — 1 and 770* = 0. The 
curved metric g^ u is then expressed in terms of a lapse function N relating 
the effective proper time to the coordinate time, a shift vector /3* expressing 
how the spatial coordinate system appears to shift among the various time 
slices, and a curved spatial metric hij with the inverse h 13 and determinant 
h. Letting g^ be the inverse curved metric as usual, we have g 00 = —N~ 2 
(the inverse metric being most convenient here), gtj = hij , and goi = hij ft 3 . 
The indices for three-dimensional quantities are raised and lowered with h l3 . 
Dropping a divergence from the Hilbert-like action above, we have the FMS 
massive version of the standard (3+l)-dimensional Lagrangian density 


c = NVh 

3 R+K ab K ob -K 2 +m 2 ( 

+m 2 





(46) 


Hereafter, we drop the superscript on 3 R. 

Defining canonical momenta as usual, we obtain the usual results 


r)C 

—— = Vh{K 13 - h l3 I<), 

ohijfi 


Pi = 


dC 


dp 


= 0, 




(47) 

The four vanishing canonical momenta are called primary constraints in the 
context of constrained dynamics m- 


10 The ADM split is a noncovariant (3+l)-representation of the metric tensor of Rie- 
mannian space first proposed by Arnowitt, Deser, and Misner. 
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Performing the generalized Legendre transformation and using the pri¬ 
mary constraints gives the canonical Hamiltonian density 


H = N 
where, as usual, 


H 0 + m 2 Vh h^rjij - 1 


+ f5 i H i -m 2 ^ + 


m 2 y/h 
2 N 


(l-r^-W), 

(48) 


Ho 


1 

y/h 



n ij 2 7F 


- VhR, 


Hi = - 2Djir 


3 
i ’ 


and Dj is the three-dimensional torsion-free covariant derivative compatible 
with hij. For m = 0, we recover the usual form that is purely a sum of 
constraints, but m/0 destroys that form and leads to six, not two, degrees 
of freedom. We note that we have retained the zeroth-order term —m 2 y/—r], 
and Minkowski space-time hence has zero energy, as it should. Boulware 
and Deser omitted this term [23]. Varying the lapse N and shift vector 
/?*, we obtain the secondary constraints, namely, the modified Hamiltonian 
constraint 


dH 


Ho + m 2 Vh 


1 + g ht3r hj 


m 2 \/h 

2N 2 


(1 - raj??) = 0 


(49) 


and the modified momentum constraint 


dH 

w 


H 2 - 


m 2 Vh 


N 


r]ijP j = 0 . 


(50) 


These constraints are second-class [48] . As Boulware and Deser pointed 
out, we can use these relations to eliminate the lapse and shift from the 
Hamiltonian density to obtain a partly on-shell Hamiltonian density purely 
in terms of the true degrees of freedom and their momenta: 


H 


12 m 2 Vh 


Ho + m 2 Vh 




+ HiHjrjd - m 2 y/^rj. 


Expressing the lapse in terms of the true degrees of freedom, we obtain 


N 2 


hm 4 

2m 2 y/h\Ho + m 2 Vh{h l i rjij / 2 — 1)] + HiHjrfi 


(51) 


Changing the variables from the lapse N to the recently popular “slicing 
density” a = N/y/h @9] allows writing the on-shell Hamiltonian density as 
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7i = m 2 (l/a — \J—r]). Boundary terms have been omitted, but they vanish 
in some cases because of the Yukawa fall-off for localized sources. 

We now consider the linearization of the exact Hamiltonian density. The 
slicing density a = N/\/h has the virtue of reducing the number of radicals 
in the off-shell Hamiltonian density and also giving the on-shell Hamiltonian 
density a simple form. Before linearization, we have 


TL = a 


7T 


tj K kl ( h ik hji - ^hijhki ) -hR + m z h\ -1 + ^-h l] ruj 


+ 


m 


+ - m 2 ^rj + — (1 - r)ij0 l f3 J ). 


(52) 


We now let rjij = 5ij , hij = 5ij + cfrij, and a = 1 + a. We are only inter¬ 
ested in the kinetic term. The quintic and quartic terms containing an 2 ^ 2 , 
n 2 cj) 2 , and all 2 (j) (indices suppressed) are dispensable, but dropping the cu¬ 
bic term o(7r^ tt 13 — /2) creates serious problems connected with the 

term — tt u 7t 33 /2 that did not arise in the exact theory. The kinetic term con¬ 
taining the worrisome wrong-sign scalar —7r 2 /2 is located in the Hamiltonian 
constraint of the exact theory, where its ability to do damage is mitigated, 
but after linearization, the —7r 2 /2 term leads to troubles, such as by radiat¬ 
ing arbitrarily much negative energy away or by permitting the radiation of 
arbitrarily much positive energy, leading to instability. If we had the cubic 
terms such as a(ir l:] ir l:] — 7r M 7r- ?J /2) and perhaps its spatial derivative ana¬ 
logue, then the approximate Hamiltonian constraint would still substantially 
resemble the exact form and might behave better. 

Boulware and Deser EU, noting the on-shell square root form of 7i above 
(and lacking the zeroth-order term), commented that “the Hamiltonian form 
(in terms of 6 degrees of freedom)... appears to be positive definite. Since 
in addition, the linearized approximation here corresponds to a scalar-ghost 
admixture, and so gives the linearized Einstein interaction in the weak-field 
limit, it would seem that this model has at least two improvements over 
[the empirically doubtful Fierz-Pauli theory and generalizations thereof]: 
Its energy is positive and it has correct linearized behavior. However, it 
is unacceptable: The vacuum is not a local minimum, but only a saddle 
point, as may be seen by considering equilibrium (static) configurations, or 
simply expanding H to quadratic order, where it is found to agree with the 
linearized (ghost) version H . That for appropriate excitations the quadratic 
part of H can be negative may seem irrelevant in view of the apparent posi¬ 
tivity of the complete H .... Unfortunately, the argument of the square root 
is not intrinsically positive... even though its positivity is required for the 
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theory to make sense, i.e., for N 2 to be positive... (otherwise, the effective 
Riemannian metric ‘seen’ by matter will become pathological). Therefore 
one would have to impose that the excitations respect this requirement, i.e., 
cut off arbitrarily those modes which take H below its vanishing vacuum 
(g = rj) state value. Instability near vacuum (g ~ rj) is the reason for re¬ 
jecting this and other models whose quadratic mass is not of Pauli-Fierz 
form.” 

This is a puzzling argument because N < 0 gives a singularity; there is 
hence no need to “impose” IV > 0 by hand. A remaining question is whether 
the theory hits IV = 0 so often that singularities form in mundane contexts. 
No argument to that effect has been given, while the fact that N —> 0 
implies gravitational time dilation suggests that N has little tendency to 
vanish [50]. There seems to be no need for Ti to be positive everywhere as 
long as H = f Ti d 3 x is positive (or nonnegative) for some suitable boundary 
conditions. The restoration of the negative zeroth-order term to Ti implies 
that Ti is bounded from below but not by zero. 

There are currently (to our knowledge) no known solutions of the non¬ 
linear field equations, exact or numerical, of FMS-Logunov or any other 
“ghost” theory of the families considered here that have negative total en¬ 
ergy. The same is true for solutions that indicate instability by radiating 
negative net energy. Given the need for a nonperturbative treatment, the 
question of stability might best be resolved with the help of numerical rela¬ 
tivists. It suffices to work in spherical symmetry, where the wrong-sign held 
can radiate but most of the right-sign fields cannot. 

8 Causality 

Given that massive gravities are considered in Minkowski space-time with 
an observable background metric rj^ u , a further issue worth considering is 
whether the null cone of the hat background metric is a bound of the effective 
curved metric g^ v . The hat metric is observable, and violation of the null 
cone of w hence implies backwards causation in some Lorentz frames, which 
is usually rejected. Causality for higher-spin theories has already caused 
trouble in the case of spin-3/2 helds. Velo and Zwanziger [51] concluded that 
the “main lesson to be drawn from our analysis is that special relativity is 
not automatically satished by writing equations that transform covariantly. 
In addition, the solutions must not propagate faster than light.” 

The argument has been made that massive gravity leads to causality 
violation in the sense of special relativity (this is relevant because of the 
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observable flat metric pi, m)- As Chugreev rightly notes, the static field 
of sources, if any, must be taken into account; for cosmological models, the 
presence of matter everywhere might suffice to preserve causality [53]. But 
surely it is a contingent rather than necessary truth that the universe is filled 
with matter everywhere f^l Gravitational radiation decays as 1/r, and the 
static field due to localized sources decays as e~ mr . Therefore, the radiation 
eventually wins and threatens the proper relation between the two null cones. 
Because we wish to regard many solutions without matter everywhere and 
with gravitational radiation as physically meaningful (although apparently 
not corresponding to the actual world), some additional strategy for ensuring 
the correct relation between the null cones of the two metrics is needed. 
The only option that comes to mind is to install artificial gauge freedom, 
perhaps using parameterization along the lines of [ 52] , [54] . and then to 
use the same strategy that we used to impose ^-causality in the massless 
case m- The resulting gauged massive gravity has a gauge transformation 
groupoid, not a group. It is noteworthy that parameterization yields results 
that for the lowest-order terms resemble Stueckelberg’s trick for introducing 
gauge freedom into massive Proca electromagnetism. Stueckelberg’s trick 
has sometimes been used in the lowest order in gravity [55], but it remained 
unclear what the generalization to nonlinear field equations might be. It 
seems plausible that other methods for installing artificial gauge freedom, 
such as the BFT procedure [56] or gauge unfixing [57], ought to give similar 
results, although we have not investigated those questions carefully. 
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